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Abstract 



The Kubo formula for the conductance of classically chaotic systems is 

analyzed semiclassically, yielding simple expressions for the mean and the 

variance of the quantum interference terms. In contrast to earlier work, here 

times longer than 0(log h^ 1 ) give the dominant contributions, i.e. the limit 

h — > is not implied. For example, the result for the weak localization 

correction to the dimensionless conductance of a chain of k classically ergodic 

1 — (fc+l) -2 , interpolating between the 

ergodic (k = 1) and the diffusive (k — » oo) limits. 
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When a micron-sized conducting sample is cooled to sub-Kelvin temperatures, the 
electrons' coherence length increases beyond the system size, and classical dynamics no 
longer applies. The conductance through the sample becomes sensitive to very small 
magnetic fields, of the order of a flux quantum through the sample. Interest in these 
well-known mesoscopic effects, known as Weak Localization [TJ (WL) and Universal 
Conductance Fluctuations || (UCF), has recently intensified, notably due to experiments on 
samples so clean that impurity scattering is no longer important ||. Originally WL and UCF 
were observed in disordered conductors with diffusive dynamics, which are described well by 
diagrammatic perturbation theory with respect to the impurity potential |l]|§. The recent 
clean systems require different theoretical tools, and calculations using Random Matrix 
Theory |4]] (RMT) and the Nonlinear a-Model || (NLoM) have appeared. These assume 
ergodic dynamics, meaning that all states at the Fermi energy are accessible "with equal 
probability". A semiclassical theory for these effects should allow for any dynamics, not 
limited to diffusive or ergodic behavior. 

Although the relevance and importance of semiclassical concepts was recognized from the 
outset ||[7| , they were used at first only for interpretation, and explicit calculations based on 
the Semiclassical Approximation (SCA) appeared only later @-[TD|. These calculations are, 
in a sense, beset by difficulties. For example, although the dependence on a magnetic field 
was obtained, the absolute magnitude of the interference effects had to be calibrated to either 
diagrammatic theory [|7| or RMT P,1Q|. The reason for this difficulty is obvious: both the 



mean (WL) and the standard deviation (UCF) of the quantum interference corrections to the 
conductance are smaller in powers of h than the leading, classical term. Such corrections to 



semiclassical propagators in chaotic systems are only very recently being studied ||11|| . In the 
perturbative theory, such corrections may be obtained by using diagrammatic elements called 
Hikami boxes ||12|| . They are known to give a vanishing contribution to the conductance, if 



it is expressed through the Kubo formula rather than the Landauer formula ||T3| (the latter 
served as the starting point for Refs. ) • 

Mesoscopic systems typically have non-separable Hamiltonians with chaotic classical 



dynamics (integrable or intermediate systems require separate consideration). The SCA for 



the Kubo conductivity in such chaotic systems has been developed by Wilkinson He was 
interested in the genuine semiclassical limit, h — > 0, and emphasized that his expressions fail 
for orbits longer than the Ehrenfest time i^, i.e. those which are long enough for an initially 
minimal wavepacket to spread due to classical chaos (Lyapunov exponents) and loose its 
correspondance to a point in classical phase-space. However, often (e.g. for the system of 
Ref. ||) the lengths of the relevant orbits, or the escape time t esc , is much longer than 
because tE — > oo only logarithmically as H — > 0. The fact that the SCA can give accurate 
results also in this "mixing" regime, was clearly demonstrated |15| only after the work of 
Wilkinson had been published. 

An application of the SCA to WL in the mixing regime will be described here. Some of 



the details, and a similar treatment of UCF are given elsewhere ||16|| . The conducting system, 
connected to two particle reservoirs through ideal leads, is described by a degenerate Fermi 
distribution of noninteracting electrons (quasiparticles) moving in an effective potential with 
completely chaotic (i.e. hyperbolic) classical dynamics. In order for the SCA to apply, 
we assume that all the features of the potential are smooth on the scale of the Fermi 
wavelength [I7|, and that the spectrum is essentially continuous, i.e. that the level spacing 
is much smaller than h/t esc . For simplicity, we describe here only the zero temperature 
DC conductance, and ignore dephasing. Spin degeneracy, time reversal symmetry, and the 
absence of additional symmetries are assumed. We concentrate on the magnitude of the WL 
correction at zero magnetic field — for the semiclassical description of the field dependence 
see Refs. @-[ig. 

The relationship j(r') = / dr a_(r, r')E(r), where j(r) and E(r) are the current density 
and the electric field in the sample, defines the space-dependent conductivity a(r, r'). The 
conductance G of a finite system is given in terms of g_(r, r') by dividing the dissipated 
power, jdr' E(r') -j(r'), by the voltage V squared: 

G = ^Jdrdr' E(r')£(r, r')E(r) . (1) 
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A straightforward derivation, starting from the Kubo formula and using the SCA, gives a 
result of the type 

£(r,r') oc J2 A a Ay^- s ^ /n ... , (2) 

a,/3e{r,r',E F } 

where a and j3 label the classical trajectories starting at r, ending at r', and propagating at 
the Fermi energy Ep. The amplitudes and actions corresponding to these orbits are denoted 
by A a and S a etc., and the dots represent further factors which depend, e.g., on the momenta 
of the orbits a and (3. The derivatives of the action S a = J a p-dr are given by = — p a , 
= p' a and = t a , where p a , p' a and t a are respectively the initial momentum, final 
momentum and duration of the corresponding orbit. 

The classical conductivity is obtained by removing all the interference terms, i.e. retaining 
only terms with a = /3. The resulting sum, J2 a \A a \ 2 ■ ■ ■, may be replaced by a mathematically 
equivalent integration over 6 functions , as follows. Any given classical initial conditions, 
r and p, and propagation time t, imply a unique final point in phase space, which we denote 
by r t and p t . The distribution of classical paths is thus f(r', p', t; r, p) = 8(r' — r t )S(p' — 
p t ). It is convenient to factor out a 5 function due to conservation of energy, defining the 
distribution on the energy hypersurface through /s(r', p', t; r, p) 5 (^H(r, p)—H(r', p')) = 
/(r', p', t; r, p) where H(r, p) is the classical Hamiltonian for the system (for an ergodic 
system f% approaches a constant at long times). In terms of this distribution, the classical 



conductivity becomes [16 
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where e is the charge of the electron, m is its (effective) mass, h is Planck's constant, d 
is the dimensionality of the system (2 for a two-dimensional electron gas, 3 for a bulk 
sample), and the momentum integrals are restricted to the Fermi surface, / dp% . . . = 
I dp 5(H(T, p)-Epj .... For a uniform system this expression can be integrated over 
r' and averaged over r, giving a more familiar form: a d = J °° dt (p(O)p(t)). Here 
v = hd 2 Yol J drdpE is the density of states (the factor of 2 accounts for the spin degeneracy) , 



and the angular brackets denote averaging over all (r, p) points on the Fermi surface [p(0) =p 
and p(t)=Pt]. 

The WL correction to the conductance is given by averaging the interference terms in 
Eq. (0). We assume that in the mixing regime S a and Sp are uncorrelated, except when 
a and f3 are related by symmetry, i.e. (3 = a T where a T is the time-reversed partner of a 
(a and a T have the same action and amplitude). The possibility of (3 = a T occurs only if 
r = r' because (3 G {r, r', E E } and a T G {r' ,r, E E }. Clearly, the ratio between the average 
conductivity and the classical conductivity increases gradually from 1 to 2 as r' approaches 
r, and we therefore consider all contributions for which (3 approaches a T when r' approaches 
r. For r' ~ r, the phase follows from the derivatives of the actions mentioned above. Using 
p'p = —p a for j3 = a T , we find that for the average conductivity a term 

O/o^ poo p 

Ag WL (r, r') ~ ^j- d J Q dtj dp E dp' E f E (r',p f , t; r, p) p' p e ^ + v')W^)/n (4) 

should be added to the classical term of Eq. (|3]) . 

The phase factor may be used to perform some of the dr and/or dp E integrations 
in the stationary phase approximation. However, the sharp 5-function nature of f E 
leads to difficulties: in the standard SCA the integrations in Eq. (f|) are performed first, 
giving back Ag WL (r, r') ~ T, a e{r,r,E F } |A*| 2 exp (i(p a +p' a )(r' -r) /ft) .... When the dr' 
integration is then performed, the phase is only stationary if a is a self-retracing orbit 
(a = a T ), and these contributions have already been included in a cl (r, r'). There are thus 
no stationary phase contributions to WL in the strict (% — > 0) sense [0. In the mixing 



regime there are many orbits with p a + p^ different from but so small that the phase 
U?Q,(r, r', E E ) — S a T(r, r', E E )^j /his negligible throughout the relevant dr' integration region, 
and it is unclear how their contribution is to be accounted for. An alternative is to replace 



f E by a smooth function f E (r', p', t; r, p) describing the density of classical orbits, which is 
obtained from the original f E by averaging over small ranges in initial and/or final conditions. 
These ranges are taken small enough (compared to the Fermi wavelength in r, and ft over 
the size of the integration region in p), so as not to affect Eq. (|4|). For times t significantly 



longer than the Ehrenfest time is a smooth function |19 . 

The integration over the components of r' and p' transverse to the direction of p can 
then easily be performed |p0 |: the stationary phase condition identifies p' with — p and the 



transverse components of r' with those of r, and the integral gives a factor of h d ~ l . The 
integration over the longitudinal component of p' gives a factor of 1/vp, where vp is the 
Fermi velocity. The phase is independent of the component of r' parallel to p (for d—2 and 
a specific orbit a, this component can in fact be defined by the condition S a = S a r). We 



therefore label the effective length of this integration region by Z(r, p), leading to our main 
result 



ml „ 2e 2 r°fj_^ tj- -^--t (p-E) 2 /(r,p) 



AG WL ~ - — dt drdp E f E (r, -p, t; r, p) -A-liZ. (5) 

h Jo J m 2 V z vf 

(the minus sign is due to the identification of p' with — p). The precise definition of 
Z(r, p) becomes apparent if one notices that the component of r parallel to p should also 
be integrated along the same segment /(r,p), and that due to the (p-EJ/mvp factors, 
each of these integrals gives the voltage difference between the two ends of the integration 
range, AV(r, p). The integration range can be curved and should extend beyond the next 
or previous scattering event in the ±p direction [F22J, but it is effectively limited by the 



"Ehrenfest length" Vptp, because of the averaging in Z(r, p) = AV(r, p)/(p-E). Note also 
that the classical results for E(r) may be used in Eq. @, and higher order corrections to 
E(r) may be ignored. 

Applications of Eq. (El) require a detailed knowledge of the classical distribution of orbits 



/g(r',p',t; r, p) and voltage drops V(r, p). For some systems this information is readily 
available — for example in a diffusive system is known from the diffusion equation, and 



l(r, p) fluctuates within the system, but is known to be equal to twice the transport mean 
free path, on the average (because both forward and backward propagation are included). 
Eq. (|5|) then reproduces the result of Ref. 0. A more interesting application is to a chain 
of k ergodic cavities, connected in series and to the two particle reservoirs through k + 1 
ideal leads of equal widths. The assumption of ergodicity within each cavity, meaning that 
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t esc is long compared to the ergodic time for that cavity, implies that /^(r', p', t; r, p) is 
independent of the fine details of the initial and final positions for long times. It also implies 
that the electrostatic potential within each cavity is a constant ||23|| , with a potential drop 



of V/(k + l) across each lead, i.e. AV(r, p) = ±V/(/c + l) for all points r within any one 
of the leads. It is convenient to define a classical dynamic probability pi, m (n), equal to the 
probability that an electron will be in the mth cavity, having started from the /th cavity 
and having traversed through a lead n times. The sum Y^=oPi,m(ji) represents the total 
number of times that an electron originating in cavity I will be found leaving cavity m, and 
is obviously related to integrals over similar to those of Eq. ([5]). Only the special case 
I — m is relevant to Eq. fl5|), because if p' = — p, and r' = r is a position within one of the 
leads (the only place where E(r) does not vanish) then the conditions (r, p) are directed 
into the same cavity that (r', p') is directed out of. In terms of Pi t m( n ), Eq. @ is 

r, 2 1 k oo 

which can easily be generalized to any configuration of ergodic cavities connected by leads 
of various widths. 

A dynamic equation, similar to the diffusion equation, may be written for pi >m {n): 

Pi, m {n+1) = - [pi >m -x{n) +pi >m+ i{ri)] , (7) 

with the initial condition pz, m (0) = 5^ m and the boundary conditions Pi${n) = and 
Pi jk+1 (n) = 0. It is easily solved: Pi, m (n) = Yh=\ ®i n A,m, where a< = cos j^ir, and 



A,i = V fc+T s ^ n fc+T 77 ' ^ s ^ n °^ ner evaluations of WL and UCF, only the eigenvalues 
— 1 < oli < 1 play a role, and the necessary summation can be done: Z)f =1 jz^: = k(k+2)/3. 
The final result for this system is thus 

AGWL = -x K 1 " jkh?) ■ (8) 

The result for UCF, which is obtained through a similar but more lengthy calculation in 



Ref. [16|, turns out to involve the 4th rather than the 2nd inverse power of k + 1. Both 



results reproduce those of RMT || for a single ergodic cavity, k = 1 (k = describes an 



ideal wire), and approach the known results for a one-dimensional diffusive wire [23| when 



k — ► oo. A similar system has been studied in Ref. |25]], using the NLcxM, and giving a set 
of more complicated results for the k = 1 to k — > oo crossover. In that study adjacent ergodic 
cavities were connected to each other through matrix elements in the Hamiltonian, rather 
than through ideal leads. 

The failure of the SCA for the Landauer formula compared to its success for the 



Kubo formula, can be traced to the fact that semiclassical evolution is not unitary, and 
the SCA for a(r, r') (like the Chambers formula for a_ d (r,r') in diffusive systems) does not 
conserve current. It is somewhat surprising that current conservation is not obeyed order 
by order in h. The issue can be clarified with a more fully understood calculation, using the 
random scattering matrices of the Circular Orthogonal Ensemble |4]] (COE). The Landauer 
formula for the dimensionless conductance of a system connected to two reservoirs through 
ideal leads with N propagating modes in each lead is g = J2iLi Sj=7v+i l^il 2 - For ergodic 
cavities, the scattering matrix Sij can be taken as a random member of the COE of 2Nx2N 
matrices. Using the unitarity of scattering matrices (current conservation), Z)f=i |5i )J -| 2 = l 
for each row or column, one may re-express sums over transmission coefficients in terms of 
sums over reflection coefficients. Use of the "actual electric field" in the sample, i.e. taking 
the potential within the cavity to be appropriately intermediate between the potential in 
the reservoirs, corresponds to making this replacement with a weight of \ for each row and 
each column and gives the "Kubo formula" : 

N 2N 

9 = 17- Y.EiE 3 \S h3 \\ (9) 

Z i,j=l 

where the "classical electric field" factors are Ei = \ for 1 < i < N and E { = -\ for 
+ 1 < % < 2N. The average conductance can be obtained from either the Landauer 
or the Kubo formula by using the known result (\Sij\ 2 ) = (1 + 5ij) / (2N + 1) (the dij term 
represents coherent backscattering) , and gives g = N 2 /(2N+1). The "semiclassical" or large 
approximation violates unitarity: (\Sij\ 2 ) ~ (l-\Sij)/2N. However, it still gives the correct 
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results to order 0(N°) if the "Kubo formula", Eq. (|^), rather than the Landauer formula is 
used: (g) = y — | + C^AT -1 ) [while only the O(N) classical term is reproduced correctly in 
the Landauer formula]. Here, as in the diagrammatic calculation, higher order corrections to 
the "propagators" or the (\Sij\ 2 ) could in principle give contributions comparable to the WL 
effect being calculated, but when their contribution is integrated with the classical electric 
field E(r), it vanishes. Unfortunately, a general proof of this is still lacking, and we are 
forced to conjecture that it holds in the SCA as well. It is expected that this "classical 
electric field" could be very useful in other calculation schemes too, e.g., that of Ref. ||. 

The appearance of new, "non-strictly semiclassical" , contributions to physical quantities 
such as the Kubo conductance in the mixing regime was not emphasized earlier, apparently 
because in most applications of the SCA the density of states was considered, and only 
periodic orbits were involved. For periodic orbits the necessity of replacing f E by f E does 
not arise, as exemplified by the UCF calculation |16[], where there are two contributions - 



one involving periodic orbits and corresponding to Wilkinson's results |T^], and the other 
appearing only in the mixing regime. It is suggested that in the mixing regime these novel 
terms may dominate the strictly h — > corrections studied in Ref. JTT] . 

The important task of demonstrating the SCA results for a given potential, rather than 
assuming a known distribution f E , is left for future research. It is emphasized that as only 
the statistical distribution of classical paths is necessary, the required numerical computation 
is not as demanding as, e.g., that of Ref. [15]. It would also be interesting to evaluate SCA 
corrections to the conductivity, e.g. due to caustics, and to contrast them with the higher 
order corrections which can be obtained diagrammatically. 
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